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Abstract 

We investigate the response of dense and hot holographic QCD (hQCD) to a static 
and baryonic electric field E using the chiral model of Sakai and Sugimoto. Strong 
fields with E > {y/XM^x) 2 free quark pairs, causing the confined vacuum and matter 
state to decay. We generalize Schwinger's QED persistence function to dense hQCD. 
At high temperature and density, Ohm's law is derived generalizing a recent result by 
Karch and O'Bannon to the chiral case. 



1 Introduction 



The AdS/CFT approach [1] provides a framework for discussing large N c gauge theories 
at strong coupling A = g\ M N c . The model suggested by Sakai and Sugimoto (SS) [2] 
offers a holographic realization of QCD (hQCD) that has Nf flavors and is chiral. For 
Nf <C N c , hQCD is a gravity dual to Nf D8-D8 branes embedded into a D4 background in 
10 dimensions. Supersymmetry is broken by the Kaluza-Klein (KK) mechanism. The SS 
model yields a holographic description of pions, vectors, axials and baryons that is in good 
agreement with experiment [2, 3, 4]. The SS model at finite temperature has been discussed 
in the deconfined phase in [5] using a black- hole metric (BH). At finite baryon density and 
in the confined phase it has been discussed in [6, 7] using a KK metric. Isospin chemical 
potential and glueball have been discussed in [9]. 

In this paper we would like to continue our investigation of the model at baryon finite 
density and temperature but in the presence of a finite baryonic electric field as recently 
discussed by Karch and O'Bannon [10] in a non-chiral model, as a prelude to understand 
transport phenomena. In section 2, we briefly outline the SS model both in the confined 
KK metric and deconfined BH metric. In section 3, the DBI action at finite baryon density 
is streamlined for both the KK and BH metrics. In section 4, we discuss Ohm's law in the 
confined or KK metric. Above a critical value of the baryon electric field E > E c the vacuum 
and the dense state are unstable against quark pair creation. In section 5, we show how this 
pair creation translates to a vacuum persistence function thereby generalizing Schwinger's 
QED result to hQCD both in the vacuum and at finite density. In section 6, we derive Ohm's 
law in the BH background, thereby extending a recent result by Karch and Bannon [10] to 
the chiral case. The vacuum instability is dwarfed by thermal pair creation in the incoherent 
statistical averaging with a treshold value for the baryonic electric field starting at zero. Our 
conclusions are in section 7. 

2 Sakai-Sugimoto model 

In this section we summarize the Sakai-Sugimoto model (D4/D8-D8 set up) for notation and 
completeness. For a thorough presentation we refer [2] for zero temperature and [5] for finite 
temperature. 

At zero temperature, the confined KK metric, dilaton 0, and the 3-form RR field C3 in 
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N c D4-branes background are given by 

{-dt 2 + 5 l3 dx*dxi + f(u)(dx*) 2 ) + ( ^\ \J±Lj + u 2 dn 



\ F A = dC s = ^e*, m = l-Vh, (1) 



where x l = x 1 ' 2 ' 3 and U(> Ukk) and ^4 are the radial coordinate and four angle variables 
in the x 5 ' 6 ' 7 ' 8 ' 9 direction. B? = 7ig s NJ 3 , with g s and / s the string coupling and length 
respectively. V4 = 87r 2 /3 is the volume of the unit S 4 and €4 is the corresponding volume 
form. To avoid a conical singularity at U — Ukk the period of St of the compactified r 
direction is set to 5r = ^-^-fp>. The Kaluza-Klein mass is 

KK 



1/2 



2tt 3CA 

'KK 



5T 2i? 3 /2 • 

The parameters R, Ukk, and g s may be expressed in terms of M KK , A(= QymNc), and Z s as 
R 3 = , U KK = l\M KK l 2 s , g s = 1 A 



2M KK 9 ™" 27rM KK iV c / s 

At finite temperature, there are two possibilities. One is the same as the zero temperature 
apart from the fact the time direction is Euclidean and compactified with a circumference 
P — 1/T. It corresponds to the confined phase which is the low temperature regime(KK). 
The other corresponds to deconfined phase of the high temperature. Its geometry contains 
the black hole and the pertinent BH background is 

3 / 2 /D\3/2 / r TT 2 



ds 2 = {f(u)dt 2 E + Sated* + (dx 4 ) 2 ) + (j?) (J?L + u 2 dnj 

e * = *(l) S/4 > F^dC^^e,, m = l-?l, (2) 
The t^-direction must be periodic 



StE = T{TF T ) = fi= T (3) 

to avoid a conical singularity. The x 4 -direction is also periodic but it has an arbitrary 
periodicity. 
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Now, consider Nf probe D8-branes in the N c D4-branes background. With U (Nf) gauge 
field A M on the D8-branes, the effective action consists of the DBI action and the Chern- 
Simons action: 



C _ qDBI i qCS 

>->D8 — >->m t o B8 , 
SE 8 BI = -T 8 J d 9 x e-* tr y/-det(g MN + 2rra'F MN ) , (4) 

where T 8 = l/((27r) 8 /g), the tension of the D8-brane, F MN = d M A N — d N A M — % [A M) A N ] 
(M,N = 0, l,-- ,8), and gMN is the induced metric on D8-branes. The expressions are 
written for the Minkowskian metric. 



3 DBI action 

The induced metric on the D8 branes from the gravity background (1) and (2) may be 
written as 

d- s m = 9ttdt 2 + g xx 8 i:j dx l dx 3 + g vu dU 2 + gssd^l (6) 



= a I - I + I - I Hijdx'dx- + 7 dU'+\-\ U'diti, (7) 

where for the KK background 

and for the BH background 

The embedding information is encoded only in 7 and thereby guu- We will use the ab- 
stract metric notations (6) to treat the confined (8) and deconfined (9) coherently in formal 
evaluation here. In the next section we will plug in the specific embedding and metric form. 
To accommodate a static baryonic electric field on D8 branes both in vacuum and matter, 
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we follow [10] to define 



A t = A t (U) , A x = -Et + h x (U) 



(10) 



With the induced metric (6) and the guage fields (10) the DBI action (4) is written as 
Sdbi = J d 4 xdU X DB i 



= -74 j dU e-+g 2 ssgxx) J \g t t\g xx guu - (2m? (fe(^) 2 + guu {A x f - M(A'J^ll) 

where 74 = (2N f)T%Vi. 2Nf comes from the fact that we consider Nf branes and anti-branes 
and V^{— 8/37T 2 ) is the volume of the unit S 4 which is due to the trivial integral over S 4 . ' 
is the derivative with respect to U and ' is the derivative with respect to t. Since (11) is 
purely kinetic, the conjugate momenta D and B are conserved. Specifically, 



_ <9£pBI _ _0 2 

= dA' ~ 6 ^ ss ^ xx 

_ (9£ DB I _ -0 2 

= dA' ~ 6 SssSxx ' 



-74{2na') 2 g xx A' t 



y/\9tt\9xx9uu ~ (27ra') 2 (9xxA' t 2 + guuE 2 - \g tt \h' 2 ) 
74(27ca') 2 \g tt \h' x 

\f\9tt\9xx9uu - (27ra') 2 (9xxA' t 2 + guuE 2 - \g u \h' 2 ) 



(12) 
(13) 



By rewriting A' t and h' x in terms of B, D and E, we have 



9xxA' t (U) 2 = 
\9tt\K(U) 2 



1 



9uu(\9tt\9xx ~ {2na') 2 E 2 ^ 



(27ra ; ) 
1 



l^\9tt\D —2 



(27ra') 



,, 2 9xxB 



74 2 (27ra'y\g tt \g xx e-^gi s + \g tt \D* - g xx B* 

9uu(\gtt\gxx - (2W) 2 E 2 ) 

7i 2 {2na'Y\g tt \gl x e-^g 4 ss + \g a \D* - g xx B* 



(14) 
(15) 



The DBI action reduces to 



S Dm = -74 I d A xdU 



/' 



p -20„4 5/21 1 1/2 1/2 

e 9ss9xx \9tt\ 9uu 



(\g tt \g xx -(2na'yE2) 



\9tt\9l x e- 



-2d> f A , Iffttl-D^ffa^-B 2 



(16) 



Notice that gtt,9xx,9ss have nothing to do with the D8 branes embedding. They carry 
information of D4 branes. Only guu carries information of the x 4 (U). It is positive for all 
U . Thus the factors outside the square root are real for all U . In contrast, the argument of 
square root may change the sign for varying U. As we will discuss below, this change in sign 
is the signal of a ground state instability or decay for large E fields. 



4 



4 Ohm's law: KK 



This decay is captured by a non-linear form of Ohm's law. For that, it is useful to change 
variable 

U = U (l + Z 2 ) 1/3 , (17) 

where U$ is the coordinate of the tip of D8-D8 branes' cigar-shaped configuration, which 
is different from C/kk in general. The range of Z is (0, oo) contrary to U whose range is 
(Uq,oo). Also this range can be extended to (—00,00) if we consider D8 branes (— 00, 0) 
together with D8 branes (0, 00) in a natural way. It enables us to deal with the ADHM 
instanton solution in R 4 [3]. It also makes the parity property of the meson fields explicit 
[2]. For completeness, we note the following useful relations 

K = l + 2?, U = U KW, dU = ^^dZ, / = l-(^) 3 ^- (18) 

From here on and for simplicity, we follow Sakai and Sugimoto [2] and choose U = Ukk- 
The DBI action then simplifies to 



>DBI 



-a / d^xdZK 1 ^, 



\l l + ^- 5/3 ' (19) 



where 



_ N c N f \*M* K 3V 

" 3%5 ' " 4A 2 M| K • [Z[)) 

In dense hQCD baryons are sourced by BPST instantons in bulk with a size of order 
They are point-like at A — > 00. Thus the DBI action and the matter sources read 

£tot = £dbi + n B 5(Z)A t (Z) + n B v x 5(Z)A x (t, Z) , (21) 

where is the baryon - anti baryon density and n B is baryon + anti baryon density. The 
first source contribution is that of static BPST instantons at Z = as initially discussed 
in [7]. The second term is their corresponding current with a velocity v x ~ l/\N c with a 
baryon mass M B ~ N C \M KK . Note that we have renormalized the field here by 1/N C 
and identified the baryon chemical potential as A^(oo) = ji B — m B [7]. 
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The equations of motion are 



D' = n B 8{Z) , B' = n B v x 5(Z) . (22) 

Thus 

1 1~ 

D = -n B sgn(Z), B = -n B v x sgn(Z) , (23) 

where sgn(Z) reflects the symmetry of D8 and D8 branes (chirality). We note that the 
conserved momenta D, B are odd functions of Z since the baryonic field is an even 
function of Z. 

For a finite baryonic electric field E, the current contribution in (21) is seen to increase 
linearly with time in the action. This is expected since the static electric field pumps energy 
in the system. For times t ~ M B ~ iV c A the present stationary (time-independent) surface 
analysis is flawed. This notwithstanding, the action variation with respect to A t yields 



$A t S t ot = I dZ 



' )L d z (SA t )+n B 5(Z)SA t (Z) 



5(d z A 



dZ (^n B sgn(Z)d z (6A t fj +n B SA t {0) 
= n B 5A t (oo) . (24) 

where we used the on-shell condition and A t (oo) = A t {— oo) = fi B — m B . Note that the 
contribution from the source term is cancelled by the boundary contribution of the DBI 
action at Z = 0. As a result the on-shell action may be considered as a functional of A t (oo) 
only and we may set A t (0) = 0. Similarly for A x (t, 0) = 0, 

£4A ot = n B v x 5A x (t, oo) . (25) 

The former is the charge, while the latter is the current. At finite density S to t plays the role 
of the grand potential. Thus 



S = -a d 4 xdZK l ' & , 



K-j^E 2 

\ 1 I »B-"B^ y-5/,s' ^ 
\ 1 + 4a 2 fe n 

on shell. For v x = E = this result is consistent with our previous result i.e. Eq.(30) in [7] 
which is indeed the grand potential. 
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For < E < E c = J x (= n B v x ) is bounded, 

J x < ^4a 2 b + n 2 B , (27) 

for S to be real. For E > E c , the numerator of (26) flips sign at 

K* = j^E 2 , Z* = ±VbE 2 - 1 . (28) 

We demand that this flip is compensated by the denominator for arbitrary v x . Using in 
the denominator we get 



- imk^ N ^{ltr Eim + nl<>iE} ' (29) 

In the unstable vacuum, the ensuing Ohmic's conductivity is 



j, i „„ r x ^ 1/3 



N c N f T7 E2/S ■ (3°) 



E 25/33^7/3 c / y Mrk 

This pair conductivity follows from quark pairs and not from baryon pairs as it scales with 
N c Nf. E c is strong enough to cause deconfinement of quark pairs. For n B ^ the second 
contribution in (29) is that of the baryons and anti baryons moving under the action of the 
strong electric field, with At> ~ Et/Ms ~ t/N c . Note that for E = 0, the minimum of (26) 
is for v x = 0. 

For E > E c both the vacuum with n B = and the dense baryonic state with n B ^ are 
unstable against pair creation of quark-antiquark states as opposed to baryon-antibaryon 
states. This is clearly seen from the threshold value E c 

M KK 2 hA-nMl . . 

E < = ^i = ^ M ™ X = ^Nf- (31) 

with M B = 8tt 2 k,M K k and k = [3] which is much smaller than M§. The baryonic electric 
field is strong enough to pair create quarks with constituent masses of order \f\M KK 1 . 



1 lt is interesting to note that in the BH background the thermal shifts of heavy quarks is tt y/XT/2 with 
T in the unconfined phase being the analogue of Mkk in the confined phase. 
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5 Persistence Probability 



The cold and dense states described by hQCD above are unstable for E > E c , meaning that 
they decay to multiparticle states that are likely time-dependent. Following Schwinger, we 
will characterize this decay through its persistence probability 



|(0 + |0_)| 2 = e- 2Im5 , 



(32) 



where ImS* is the imaginary part of the action S (26). For finite n# and v x = 0, the action 
S reads 



S = -a J d 4 xdZ(l + Z 2 ) 1/6 



Z 2 + 1 - £ 5 



Z 2 + l 



1 + + Z 2 )- 5 / 3 V 1 + jV \\ + z 2 )- 5 / 3 



(33) 



with £ z 



b T?2 



Ml 



E , jV = and after regularizing the action by subtracting the £ = 
contribution. E c corresponds to £ c = 1. For £ < 1 the action S is always real, but for 
£ > 1 the action develops an imaginary part from the integration interval (— Z c , Z c ), where 
Z c = Vfi 2 - 1. Thus 



ImS = ±a I d A x dZ(l + Z 2 y/\ *L_^(£ _ i) . (34) 



For — the integrals unwind analytically 



ImS* = ±air d x 



I 



(£ 2 -!),/•', ( -|,|,2,1 - £ 2 ) 1) 



(35) 



where 2 Fi is the hypergeometric function and has the asymptotic behaviour as follows. 



.,/•, ( --,-,2,1 - £ 2 
' 6 2 



i + l(£-i) + I ^(£-i) 2 + 



r(2/3) 

0Fr(13/6) 



£1/3 + 



2r(-2/3) 1 
v^Fr(-l/6)£ 



+ 



(£-1) 
(£ » 1) 



(36) 
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The persistence function is then 



|(0+|0_)| 2 = e ~ a ' (£2_1) 2fi(-ii2,i- £ 2 )e( £ -i)^ 

= 1 (£<1) 

e -a'[2( £ -l)+1.17( £ -l) 2 +-] ( £ ^ 1 \ 

e -a'[0.71 S+-] (£ ^ ^ (37) 

with a' = got J d 4 x = NcN ^ w ^ KK J d 4 x, after chosing the negative sign for decay. 



6 Ohm's law: BH 

Since the vacuum decay under large E"s so does the coherent finite baryonic state. But 
what about the finite temperature problem? As finite temperature involves a statistical 
ensemble averaging, we may suggest that the unstable ground state is statistically irrelevant 
and proceed to analyse the effects of a baryonic field on the excited states (unstable by fiat) 
in the ensemble average. This will be checked a posteriori below. 

In the BH background there are two possible gravitational configurations: 1/ a U-shaped 
(chirally broken phase) and 2/ a parallell-shape (chirally symmetric phase). The former 
yields U bounded from below by Uq. The combination guQxx has a positive minimum so the 
numerator is always positive for sufficiently small E. The nature of the transition which is 
suggestive of a metal-insulator transition [11] will be discussed elsewhere. 

For high enough temperature the stable configuration is not the U-shaped configuration 
but the parallel configuration which is connected to the black hole. i.e. ^ = 0. Our intial 
instanton sources have now drowned into the BH horizon. So the ensuing analysis is the 
same as in the D3/D7 model [10], with the general formula of the conductivity for Dq/Dp 
given ((5.7) in [10]). Here and for completeness, we compute the conductivity for the parallel 
D8-D8 branes set up in the BH background. 

We only need to consider the positivity condition for the argument of square root as 
before. As U — > Ut both the numerator and denominator are negative since gu — > 0. As 
U — > oo both the numerator and denominator are positive. So by choosing B,D,E we can 
choose the numerator and denominator in (16) to flip sign for the same value U = U* [10]. 
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For the numerator 



u=u, 



\9tt\9xx 

U* = (C/| + fl 3 (2ro'E) 2 ) 1/3 • (38) 



Inserting this value of U* in the denominator yields the induced current 



Jl = (^a'f\9u\9lxe'^9 4 ss + ^J 2 t 

9xx 



u=u. 



where J x = B and J t = D (= rig) are now defined as in [10]. Setting U T = ^-T 2 R 3 , 
A = g s N c yield the Ohmic conductivity for the chiral SS model 



where 



J x j f 4l s NfN c \T 2 ^ 2 d 2 

a= E = n-^27 — ) (i+e)/ + tt^' (40) 



d = ^, » ( 41 ) 



2 5 tt 3 T 3 A/ s ' 2 5 tt 3 T 3 A/ s 

which is consistent with the result in [10] for massless but non-chiral quarks. The induced 
thermal current sets in for any E > (large or small) with a conductivity a of order 
N c Nf\T 2 l s at high temperature, which involves only thermal pairs with zero treshold for 
E. It dwarfs the induced vacuum pairs by a factor of A 2 / 3 . The unstable vacuum state is 
statistically irrelevant. This is not the case at T = and/or very large baryonic densities. 



7 Conclusions 

We have extended our recent holographic analysis of the SS model at finite density, to the 
case of finite temperature and finite baryonic electric field. For E > E c the stationary SS 
ground state breaks down by quark pair creation. This phenomenon permeates both the 
cold and hot states of hQCD. The vacuum persistence probability is derived, generalizing 
Schwinger's QED result to hQCD. At finite temperature, the baryonic electric field yields 
a thermal conductivity at finite temperature and density that is a direct generalization of 
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Karch and Bannon's Ohm's law in the chiral model. We have argued that the vacuum 
instability is statistically irrelevant in hot hQCD. 

8 Acknowledgments 

The work of KYK and IZ was supported in part by US-DOE grants DE-FG02-88ER40388 
and DE-FG03-97ER4014. The work of SJS was supported by KOSEF Grant R01-2007-000- 
10214-0 and the SRC Program of the KOSEF through the CQUEST with grant number 
Rll-2005-021. 

9 Note added 

While typing our results, a recent analysis appeared in [12] that addresses similar issues in 
the model with cusp surfaces at L ^ ir. Our results are all for cusp less surfaces with L = n 
for the KK background, and the parallel or deconfined configuration for the BH background. 
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